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We consider combining two important methods for constructing quasi-equilibrium initial data for 
binary black holes: the conformal thin-sandwich formalism and the puncture method. The former 
seeks to enforce stationarity in the conformal three-metric and the latter attempts to avoid internal 
boundaries, like minimal surfaces or apparent horizons. We show that these two methods make 
partially conflicting requirements on the boundary conditions that determine the time slices. In 
particular, it does not seem possible to construct slices that are quasi-stationary and avoid physical 
singularities and simultaneously are connected by an everywhere positive lapse function, a condition 
which must obtain if internal boundaries are to be avoided. Some relaxation of these conflicting 
requirements may yield a soluble system, but some of the advantages that were sought in combining 
these approaches will be lost. 
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I. INTRODUCTION 

Binary black holes are expected to produce some of 
the strongest gravitational wave signals that lie within 
the pass-bands of both the new generation of terrestrial 
gravitational wave observatories, such as LIGO, GEO, 
TAMA, and VIRGO, and the space-based LISA detec- 
tor. All aspects of the orbital decay and signal history of 
such decaying binaries are of interest, both for enhanc- 
ing the prospects of source detection and for determining 
source parameters 0, . With this need in mind, an in- 
tense worldwide theoretical effort is being made to model 
different phases of the merger of binary black holes. 

Considerable interest attaches to the phase between 
early inspiral and loss of orbital stability and plunge. 
For astrophysical reasons many binary black holes are 
expected to enter this phase in quasi-circular orbital 
motion, with the residual orbital eccentricity following 
merely from the gradual radiative loss of energy and an- 
gular momentum. Hence, early on these quasi-circular 
orbits are well approximated by positing circular motion 
and ignoring radiation reaction. 

The first numerical models of binary black holes in 
quasi-circular orbit were obtained by Cook 0] (hereafter 
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C94), who used an effective-potential method to identify 
those particular data sets that most closely approximated 
binaries in circular orbits. The initial data were con- 
structed using the conformal-imaging formalism. In that 
approach, the momentum constraints decouple from the 
Hamiltonian constraint and can be solved analytically, 
yielding the Bowen-York conformal extrinsic curvature 
4J. The conformal-imaging approach takes as given a 
conformally-flat metric, a maximal time slice, and a two- 
sheeted topology where two identical universes are con- 
nected by one or more black holes. The identification of 
the universes manifests as an isometry condition that re- 
stricts the form of the solutions. This isometry condition 
can be used to obtain boundary conditions on the black- 
hole throats. In particular, a minimal surface boundary 
condition is obtained for the conformal factor and the 
Hamiltonian constraint need only be solved on one of 
the two isometric sheets of the hypersurface. The solu- 
tion of the Hamiltonian constraint is (near completely) 
restricted to the exterior of the black holes in one of the 
two universes and the calculation avoids both coordinate 
and physical singularities in the black hole interiors. 

Brandt and Briigmann [B| suggested an alternative 
puncture method for the construction of spacetimes con- 
taining multiple black holes. The momentum constraint 
is again solved analytically on the conformal manifold, 
yielding the Bowen-York extrinsic curvature. The confor- 
mal factor, however, is split into a sum of an analytically- 
known singular term and a correction. The analytic 
term corresponds to the conformal factor of superposed 
static black holes in isotropic coordinates and captures 
the topological character of the black holes. The Hamil- 
tonian constraint is then written as an equation for the 
correction term, which exists for non-vanishing extrin- 
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sic curvature. In this method, though, isometry across 
the throats is not assumed or used. Instead the compu- 
tational domain is all of Ft 3 minus the puncture points 
where the conformal factor is singular. However, by con- 
struction, the terms in the Hamiltonian constraint are 
regular everywhere, allowing the punctures to be ignored. 
By eliminating the need to excise the black hole interi- 
ors and to impose boundary conditions on the black hole 
throats, the puncture method provides a much simpler 
computational domain. Baumgarte Q (hereafter BOO) 
adopted this approach to construct binary black holes 
in quasi-circular orbits and found very similar results to 
those of C94. This is not surprising, as the resulting 
spacetimes only differ in terms of the underlying topol- 
ogy (the puncture method gives rise to N + 1 asymptoti- 
cally flat regions with N throats), and the effects can be 
estimated to be fairly small. 

Gourgoulhon, Grandclement and Bonazolla 0, Q 
(hereafter GGB02) adopted a conformal thin-sandwich 
decomposition of the initial value problem (e.g. ITfij) 
to construct binary black hole data. Unlike the methods 
described above, the conformal thin-sandwich decompo- 
sition has a link to dynamics. This is seen in the appear- 
ance of the lapse function and shift vector in the con- 
formal thin-sandwich equations, and in the fact that the 
time rate of change of the conformal three-geometry is 
part of the freely-specifiable data. In the conformal thin- 
sandwich approach, the momentum constraint does not 
yield the Bowen-York extrinsic curvature as a solution. 
Instead, it becomes an elliptic system on the shift vec- 
tor, which once solved indirectly determines the extrinsic 
curvature. As before, the Hamiltonian constraint is an 
elliptic equation for the conformal factor but the confor- 
mal thin-sandwich approach must be supplemented with 
a condition on the lapse. GGB02 adopted maximal slic- 
ing, which is a natural choice for the construction of equi- 
librium data. They also take the time rate of change of 
the conformal three geometry to vanish which again is a 
natural choice for equilibrium data. GGB02 also follow 
C94 in assuming an isometry across the throats, which 
allows the black hole interiors to be excised. While this 
approach is more involved than those of C94 and BOO, 
it has been suggested that the conformal thin-sandwich 
decomposition is more natural for the construction of bi- 
naries in quasi-circular orbit (e.g. 01)' an d m fact the re- 
sults of GGB02 agree more closely with post-Newtonian 
calculations than do those of C94 and BOO |llill2l ]. There 
is furthermore good evidence that the differences between 
the results of GGB02 and those of C94 and BOO are 
due in fact to the differences in initial value decompo- 
sitions, inparticular the behavior of the extrinsic curva- 
ture Ha Hi Hi - 

As we will discuss in more detail in Section IIIII the 
imposition of the isometry on the throats leads to a tech- 
nical problem in the conformal thin-sandwich approach. 
GGB02 had to regularize their solutions on the throats in 
a way that introduced a small inconsistency (i.e. a small 
violation of the constraints) into their solutions. It is 



thought that this problem may be circumvented with a 
different set of boundary conditions on the throats ^3 . 
This issue, however, suggests that one might try to com- 
bine the thin-sandwich approach with the puncture ap- 
proach, with the goal of eliminating the need for interior 
boundary conditions. 

We have attempted to combine the conformal thin- 
sandwich and puncture methods, but have concluded 
that this approach leads to some conceptional problems 
and some likely serious numerical obstacles. As we will 
argue below the approach requires several simultaneous 
conditions: (1) The construction of quasi-equilibrium 
data requires the existence of a quasi-stationary coordi- 
nate system. (2) The puncture method requires phys- 
ical fields in the vicinity of each puncture such that 
each puncture represents a separate spatial infinity on 
a multiply-connected manifold (i.e., the slice(s) must 
thread each wormhole and avoid the physical singular- 
ities). (3) Finally, the conformal thin-sandwich decom- 
position (usually) requires that the lapse be everywhere 
positive, to avoid the existence of internal boundaries ne- 
cessitating regularity conditions. Unfortunately, even in 
Schwarzschild spacetime such a slicing cannot exist, sug- 
gesting that the construction of black hole equilibrium 
data with a combined conformal thin-sandwich and punc- 
ture method is either impossible or will involve numerical 
complications that it was initially intended to avoid. 

The remainder of the paper is organized as follows. In 
Section El we review briefly the conformal thin-sandwich 
decomposition. Section IIIII shows the choice of freely- 
specifiable data that represents quasi-equilibrium, giving 
the form of the conformal thin-sandwich equations that 
must be solved. Section llVl then briefly reviews the orig- 
inal puncture method. Based on that discussion, Section 
presents how we argue that singular terms can be fac- 
tored out of the conformal factor and lapse function in the 
conformal thin-sandwich equations of Section IIIII allow- 
ing those equations to be solved as if the manifold were 
free of punctures. Freedom exists in the choice of lapse 
and the remainder of this section focuses on the poten- 
tial inconsistencies or numerical difficulties in solving the 
momentum constraints. Finally, Section IVII summarizes 
our conclusions and suggests some avenues for further 
exploration. 

II. THE CONFORMAL THIN-SANDWICH 
DECOMPOSITION 

The spacetime line element can be written in the 3+1 
form 

ds 2 = -a 2 dt 2 + -f ij (dx i + f3*dt){dx j + (3 j dt), (1) 

where a is the lapse function, (3 l the shift vector, and 
jij the induced metric on a spatial slice S. The spatial 
metric Yy and the slice's extrinsic curvature Kij satisfy 
the Hamiltonian constraint 

R + K 2 - KijK ij =0 (2) 
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and the momentum constraint 



V,/\'° - -''V,/\ = 



(3) 



on each slice E. Here Vj, Rij, and R = J tJ Rij are the 
spatial covariant derivative, Ricci tensor, and scalar cur- 
vature associated with 7^ (we have no occasion in this 
paper to use the analogous four-dimensional quantities). 
We have assumed vanishing matter sources {T^ v = 0). 
The metric and extrinsic curvature evolve in accordance 
with 



Vi/3j + Vjfa 



and 



d t Ki 



a{R l3 - 2K tk K) + KKij) 
+ K lk V ] p k + K k3 V l p k . 



(4) 



(5) 



The construction of initial-data solutions requires spec- 
ifying the spatial metric and extrinsic curvature that 
satisfy the constraint equations J2J) and (0. The ap- 
proach based on conformal decomposition |l6( provides 
a straightforward process to solving these coupled equa- 
tions. The spatial metric is conformally scaled by 



(6) 

(see 03 EH). Here^ is the conformal factor and %j is the 
conformal or background metric. The extrinsic curvature 
Kij is split into its trace K and the traceless part Ay, 
and the latter is then also conformally transformed 

K i:j = A*, + -jijK = ip~ 2 A i3 + -ip%jK. (7) 

There is no conformal scaling of K . The Hamiltonian 
constraint @ can now be written as 



V 2 ip 



-ipR 
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where V, E44, and R — 7^ Ri 



^j- 7 A t3 A^, (8) 
denote the covariant 



derivative, conformal Ricci tensor, and curvature scalar 
compatible with the conformal metric 7^ . 

In the conformal thin-sandwich decomposition, both 
the background metric 7y , its time derivative 



'tHj, 



(9) 



and the trace of extrinsic curvature K are considered 
freely-specifiable quantities. The traceless part of the 
evolution equation Q for 7^ then yields an equation for 



(L/3)i, - iiij 



(10) 



where L is a longitudinal operator whose action (L/?)* 3 
on the shift yields the symmetrized trace-free gradient, 



(L/3) ij ' = V*jS* + V j f3 l - ^f j V k /3 h 



(11) 



and where we have introduced the densitized lapse (or 
slicing function) 

q = >ip- 6 a. (12) 
Inserting (|1C)|> into the momentum constraint yields 



4 . 
3 ( 



A^-CL^Vjlna = ^atp^K+aVj [ ^u t] ) , (13) 



which can be regarded as an elliptic system for the shift 
vector f3 l . Here the operator Al (a vector Laplacian) is 
the divergence of the longitudinal operator L, with the 
effect 

A L /3 l = VjQLpf = V 2 /3 J + \v\\J k p k ) + mp. (14) 

3 J 

The shift vector components are conformally invariant 
(/?' = /3 1 ). Equations JHJ) and (|13|l are incomplete in that 
we have yet to specify how the lapse a is chosen. The 
condition on the lapse is intertwined with the choice of 
free data that yield a quasi-equilibrium state, and we take 
up these issues next. 



III. QUASI-EQUILIBRIUM DATA IN THE 
THIN-SANDWICH APPROACH 



The construction of binary black holes in quasi-circular 
orbit assumes that the emission of gravitational radiation 
is negligible, so that the spacetime is in quasi-equilibrium 
when viewed from a co-rotating reference frame. The 
time coordinate of this reference frame is an approximate 
helical Killing vector |l9|. 

Our goal is therefore to construct initial data which, 
when evolved in such a co-rotating reference frame, 
would lead to metric components that are independent of 
time. The conformal thin-sandwich decomposition pro- 
vides a natural framework for the construction of quasi- 
equilibrium data, since we can explicitly set the time 
derivative of the conformal metric to zero 



Uij = 0. 



(15) 



We follow C94, BOO and GGB02 and assume conformal 
flatness, i.e. 



li'j ~ ft] ' 



(16) 



where /y is a flat metric. It is furthermore reasonable to 
assume that the time-derivative of the trace of the extrin- 
sic curvature vanishes, dtK — 0. With this assumption, 
the trace of (0 combined with J5J yields 

V 2 (^ 7 5) = a^ 7 Q^" 8 A y A^ + ^ 4 K 2 ^j + 

(17) 

In this equation, A lJ is now given by 

A ij = -l(L/3) y '. (18) 
2a 
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We could further follow GGB02 and choose maximal slic- 
ing, K = 0, but the argument is more general and we 
merely assume that a sufficiently smooth choice for K 
is made that is consistent with asymptotic flatness and 
spatial slices. 

Our choices of 7y = fy, Uij — and K completely 
determine all freely-specifiable variables. With these 
choices, the Hamiltonian constraint JSJ reduces to 

- --ii-'AijAv + ^ 5 K 2 (19) 
and the momentum constraint i|13|) reduces to 



4(L/3) y ' 
a 



A L /3 ,; - (L/3) ij 'Vjln5 



3 v 

(20) 

Equations (|17|l through l|20[l together with appropriate 
boundary conditions now determine the conformal factor 
ifi, the lapse a and the shift l . This set of equations was 
derived independently by |2fj l2l| and has been used in 
several applications, in particular for the construction of 
binary neutron stars (see [loll^ for recent reviews). 

The choice of u%j — and dtK — guarantees a lim- 
ited sense of stationarity. However, there is no assur- 
ance that the other components of the extrinsic curva- 
ture or the conformal factor have a vanishing first time 
derivative. (Were this the case, it would imply exact 
equilibrium.) Furthermore, there is no guarantee that 
these time derivatives are even small, in some appropri- 
ate norm. There are a couple of issues to consider. 

Firstly, the boundary conditions on the problem must 
be made consistent with the limited stationarity condi- 
tions given above. Boundary conditions can be imposed 
in different ways (e.g., in the puncture method via free 
parameters at the punctures-see below) but must be re- 
flective of two black holes in instantaneous transverse 
motion or vanishing motion in a rotating frame. Even 
then, however, parameter freedom exists and only some 
parameter choices will lead to black holes in instanta- 
neous force balance as seen in the rotating frame. Most 
parameter choices will imply an acceleration of the holes 
away from or toward each other, which will be reflected 
in the time derivative of Aij . All methods to date involve 
a parameter search and some equilibrium criterion, such 
as an effective potential minimization or a virial theorem 
argument. 

Secondly, given a solution for tp, a and (3 l following 
such a parameter search, one could a posteriori compute 
the right hand side of J5J and the trace of J3) to check 
whether these time derivatives are sufficiently small in an 
appropriate sense. If they are small, then one has in fact 
constructed a quasi-equilibrium orbit. If, however, these 
time derivatives are not small, the initial data may still 
represent a quasi-equilibrium orbit but in a dynamical 
slicing. Even on a stationary spacetime with a timelike 
Killing vector, an arbitrary choice of the lapse or shift 
would lead generically to apparent time-dependence of 
the metric. In such a dynamical slicing the Killing vector 



representing the symmetry (a helical Killing vector in the 
case of binaries) would have an unknown functional form 
and it is unclear how the symmetry would be established 
or exploited. For the construction of quasi-equilibrium 
orbits it therefore seems crucial that the solution be rep- 
resented in a stationary slicing. The slicing condition 
(|17|l may be consistent with stationarity but stationar- 
ity also depends on boundary conditions. It is also an 
open question what advantage might arise in specifying 
some K ^ but any such choice would also need to be 
consistent with stationarity. 

In constructing their binary black hole solutions, 
GGB02 employed a method that is essentially the con- 
formal thin-sandwich decomposition as outlined in Sec- 
tion ^ and this section. They demanded additionally 
that their solutions be inversion symmetric, i.e., that the 
isometry conditions across the throats be satisfied as in 
C94. One choice for inversion symmetry of the lapse is to 
require, as did GGB02, that a = on the throats. This 
condition clearly leads to problems in constructing the 
extrinsic curvature via Eq. i|18[l unless (L/3) y = on the 
throats as well. This latter condition is not the natural 
isometry condition on the shift, nor is it the boundary 
condition used in GGB02. Rather, GGB02 enforced this 
extra condition, over-determining the shift, and leading 
to a small inconsistency in their numerical solution. It 
may be that these problems can be resolved by adopting 
an alternative set of boundary conditions representing 
black holes in equilibrium |l5j . 

Alternatively, one might hope to avoid interior bound- 
ary conditions altogether by combining the conformal 
thin-sandwich approach with an extension of the original 
puncture method. To set the stage for that discussion, 
we first review briefly the puncture method. 



IV. THE ORIGINAL PUNCTURE METHOD 

The puncture method can be motivated by considering 
the Brill-Lindquist 23] solution for multiple black holes 
at a moment of time symmetry (K — 0,Aij = 0). In 
isotropic coordinates {x 1 }, the conformal factor can be 
written as 



1> 



N 

E 



Mr 



2x- 



N 

E 



Mr 

2r n 



(21) 



This is a solution of the Hamiltonian constraint 1|19|) for 
TV black holes at locations x„, each parameterized by a 
mass M n . For convenience we define 



1 N li A 

1 = Mn 
- ~ ^ 2r n 



(22) 



so that the above solution is simply ip = 1 + l/p. Brandt 
and Briigmann |5( noted that for black holes with non- 
zero linear and angular momentum (A^ ^ 0), the con- 
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formal factor can be written as 



V> = H \-u, 

P 



(23) 



where u is a correction term (our definition of u corre- 
sponds to u — 1 in [f| and BOO) . Having factored out 1/p, 
the Hamiltonian constraint becomes an equation for u, 



V 2 u 



-^- 7 A ij A i ^ + ^ b K 2 . (24) 



Provided the trace-free extrinsic curvature A 13 diverges 
at each puncture no faster than 1/r 3 , the term ip~ 7 Aij A 13 
will be everywhere finite. A similar restriction on the 
second term on the right hand side of (|24(l (Brandt and 
Briigmann had assumed K = 0) requires that K van- 
ish at least as fast as r 3 at each puncture. For such 
choices of A 13 and K, u will be at least a C 2 function 
everywhere. This allows the Hamiltonian constraint to 
be solved as if the conformal manifold were H 3 , ignoring 
the punctures. In this way, there is no need to provide 
inner boundary conditions on minimal surfaces (throats) 
or apparent horizons 0, ^| . This greatly simplifies the 
solution procedure. BOO applied the puncture method 
with a (modified) Bowen-York extrinsic curvature (the 
terms that guarantee isometry are dropped for the 
reason given above). Baumgarte's initial value problem 
consisted of this analytic solution of the momentum con- 
straints and a numerical solution of Thus he was 
able to redo the calculation of C94 with the isometry as- 
sumption replaced by the puncture approach, construct- 
ing a sequence of black-hole binaries in quasi-circular or- 
bit. Quite similar physical results were obtained. 

The puncture approach relies upon the spatial slice 
E passing smoothly through each wormhole and reach- 
ing a separate spatial infinity associated with each black 
hole. The particular form of the unbounded growth in 
tp (assuming u is finite) as each puncture is approached 
x — > x n is indicative of separate compactified asymptotic 
infinities. The intent is to avoid physical singularities, 
which could be verified by computing a Riemann cur- 
vature scalar or, as is typical, using a set of coordinate 
transformations (inversion through spheres about each 
puncture) and verifying that ip and Ay are consistent 
with asymptotic flatness at the punctures. 

With the success of the puncture method in yield- 
ing results similar to the conformal imaging approach, 
it seemed worth considering a combination of the punc- 
ture method and the conformal thin-sandwich formalism. 
Such a method has the potential of conferring computa- 
tional simplicity and yet avoiding the ambiguities in the 
boundary conditions on the throats in the calculations 
of GGB02. The next section outlines the idea and its 
potential shortcomings. 



V. A COMBINATION OF THE CONFORMAL 
THIN-SANDWICH FORMALISM AND THE 
PUNCTURE METHOD 

In order to apply the puncture method to the confor- 
mal thin-sandwich decomposition, as a first step we need 
to examine how the lapse equation 117|) and lapse func- 
tion a behave in the vicinity of the punctures. From 
the similarities between the lapse equation (|T7jl and the 
Hamiltonian constraint |(SJ| we can anticipate splitting 
ai/j 7 into a sum 



2r„ 



(25) 



where the C n are constants and v is a correction. The 
term v is expected to be finite and sufficiently smooth 
that the domain once again can be taken as Euclidean. 
As with the Hamiltonian constraint, for this to be true 
conditions must be imposed at the punctures on the be- 
havior of A 13 , K, and now /3\ 

The possible values of the C n correspond to different 
boundary conditions on the lapse at the other N asymp- 
totic infinities. Each unique set {C n } gives rise to a dif- 
ferent development between time slices. The full general- 
ity of these boundary conditions need not be considered 
here. It suffices for our present purpose to distinguish 
two principal classes. 

To understand these two choices, consider a single 
Schwarzschild black hole described by isotropic spatial 
coordinates. One might choose to describe the spacetime 
with static Schwarzschild time slices, for which the lapse 
function is 



1 - M/2r 
1 + M/2r ' 



(26) 



Given that -0 = 1 + M/2r in these coordinates, we have 



M 
2? 



(27) 



which satisfies equation (|17|l and the lapse serves to link 
successive maximal slices. As is well known, these slices 
are frozen on the throat (r = M/2), have negative lapse 
on the second sheet, and a — > — 1 at the other spatial 
infinity. In analogy with the arguments of Section IIVI 
this suggests that the generalization 



atp 7 = 1 — — + v 
P 



(28) 



be used for multiple black holes. In this case the lapse 
would also be negative on each of the other N sheets and 
limit on —1. 

Alternatively, one might choose to use the opposite 
sign on the singular terms, to wit 



1 



1 



(29) 
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For a single Schwarzschild hole this reduces to a-0 7 = 
1 + M/2r and hence to a = 1. The unit lapse does 
connect two slices in the foliation of maximal slices of 
Schwarzschild found by Estabrook et. al. (24) . Each pair 
of slices in that foliation has a > everywhere, though 
in general 9,« / 0. In applying this boundary condition 
to two or more black holes, we certainly expect did ^ 
but, with a sufficient bound on v, will have a > and 
a — > 1 at infinity on the other sheets. 

These two choices illustrate the primary dichotomy in 
the choice of boundary conditions on the lapse: Either 
the lapse reverses sign in the interior of the black holes, 
which in the puncture method means that each puncture 
is surrounded by a two-surface on which a = 0, or the 
lapse remains positive everywhere, including on each then 
dynamic throat. In a nutshell, therein lie the potential 
shortcomings of combining the puncture method with the 
conformal thin-sandwich approach. 

In the former case Q28fl. the vanishing of the lapse on 
some surface around each puncture is very problematic 
for the solution of equation (|2(J[1 or for the construction 
of the extrinsic curvature from equation JTSJ. Examina- 
tion of H20|l reveals that a = is a critical surface for 
the differential equations. A regular solution for (3 l may 
be obtained by imposing the regularity condition that 
TOj(L/3) y = on this surface, where rrij is the normal 
to the level surface a = 0. However, even if a regular 
solution for (3 l can be obtained, handling such an in- 
ternal boundary is precisely the type of task that the 
puncture method seeks to avoid. Furthermore, a regular 
solution for (3 l does not ensure a finite extrinsic curva- 
ture, as reference to equation (|18|l indicates. Stipulating 
the vanishing of all five algebraically-independent compo- 
nents of (L/?) 1 -? on the critical surface appears to entail 
more freedom than is available. This is the difficulty that 
GGB02 faced in requiring both isometry on the shift and 
regularity in A lJ . 

So, if the use of the boundary conditions implied by 
equation (|28|l is doubtful, what of the latter choice ex- 
pressed by equation (|2"5|l ? Here the problem is more sub- 
tle. As we discussed in Section UTTl it seems undesirable 
to construct a quasi-equilibrium solution through some 
means that leads to a dynamical slicing. At minimum, 
this would make it difficult to verify whether or not one 
has in fact constructed a quasi-equilibrium solution. This 
requisite and the previous one lead to the question of 
whether it is possible to find a pair of slices (maximal 
or non-maximal) on which the fields appear stationary, 
that give rise to a lapse that is positive everywhere, and 
which end at spatial infinity on all sheets. The second 
and third conditions are requirements on the applicability 
and usefulness of the puncture method. The first condi- 
tion attempts to ensure that the quasi-equilibrium nature 
of data constructed with the thin-sandwich formalism ac- 
tually appear in equilibrium. As we will argue below, 
such a slicing does not even exist for a Schwarzschild 
spacetime. 

The argument is easy to make. Consider a t equal 



constant surface of Schwarzschild and a = 1 consistent 
with (|29|l . A quick evaluation of the evolution equation 
© indicates that the time derivative of A y will be non- 
vanishing. In crossing the hole with a positive lapse we 
necessarily begin to see the pinch-off of the throat. 

This argument can be made more general by consider- 
ing a standard Penrose diagram of Schwarzschild, as in 
Figure Q The curved lines with arrows represent the 
Killing flow associated with the Killing vector that is 
timelike in the black hole exterior. The Killing vectors 
are tangent to the r = const curves. The two dashed 
lines represent two spacelike hypersurfaces that have no 
special properties other than ending at spatial infinity in 
both universes. In particular, there is no reason that they 
should be maximal. In this example, the upper surface 
represents a slice that stems from evolution off the lower 
slice with a lapse that is everywhere positive. 




FIG. 1: Penrose diagram of the Schwarzschild spacetime. The 
curved solid lines are a Killing flow and the dashed lines are 
spacelike hypersurfaces. 

Figure [21 focuses on the black hole interior and just 
one of the hypersurfaces. Here we begin with no a priori 
assumption about the second slice or whether the lapse 
changes sign or not in crossing the black hole interior. 
Point "B" represents the minimal surface or throat in the 
hypersurface, as it corresponds to the smallest value of 
the radial coordinate (e.g., areal coordinate). Any slicing 
that ends at spatial infinity in both universes must have 
such a point. 




FIG. 2: Penrose diagram of the black hole interior. Point B 
represents the minimal two-surface. 

This leads to the following crucial observation. As we 
are interested in stationary slicings, the time four- vector 
t M (d/dt) for evolving to the next slice should be coinci- 
dent with the Killing vector. The spatial slice is tangent 
to the r — const surface at the throat. But the Killing 
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vector is tangent to the r = const surface everywhere. So, 
in particular, the Killing vector is tangent to the spatial 
slice at the throat. The time four-vector can be written 
in terms of the lapse and shift as i M = an 11 + /3 M . At 
the location in question, since the Killing vector is tan- 
gent to the hypersurface, the lapse must vanish at the 
throat. This is unavoidable and establishes that we can- 
not choose a slicing that connects two spatial infinities, is 
time independent, and has an everywhere positive lapse. 
The lapse must at least vanish at the minimal surface if 
the slicing is stationary. 

The argument can be extended to show that if station- 
arity is assumed the lapse must become negative interior 
to the throat. At point "C" in the figure the Killing vec- 
tor points in the upper right direction so the lapse must 
be positive there. But, at point "A", the Killing vector 
points in the lower right direction, so the lapse will be 
negative at "A" . Alternatively, if the lapse is positive def- 
inite, the slice must be non-stationary at least in some 
neighborhood of the black hole throat. 

An added consequence of these arguments is that a slic- 
ing of Schwarzschild that is stationary and has an every- 
where positive lapse can never be tangent to an r = const 
surface. This requires that one end of the slice terminate 
at the physical singularity, which obviates the application 
of the puncture method. 

VI. SUMMARY 

Several arguments point to the conformal thin- 
sandwich decomposition as a promising approach for con- 
structing binary black holes in quasi-circular orbits (cf. 
GGB02). However, assuming an isometry and using the 
black hole throats as inner boundaries leads to some 
mathematical inconsistencies (cf. Ref. ^H). This expe- 
rience suggested that one might combine the conformal 
thin-sandwich approach with the puncture method 
which seeks to avoid interior boundaries altogether. 

Combining the conformal thin-sandwich decomposi- 
tion with the puncture method for the construction of 
equilibrium data imposes three requirements on the slic- 
ing. Equilibrium data can most easily be identified when 
they are presented in stationary slices. The puncture 



method also requires that the slices connect the spa- 
tial infinities of the multiple black holes, as opposed to 
ending on a physical singularity. Finally, the conformal 
thin-sandwich decomposition appears to require an ev- 
erywhere positive lapse, to ensure a finite extrinsic cur- 
vature and, at minimum, to avoid internal boundaries. 
However, even a Schwarzschild spacetime admits no sta- 
tionary slicing that has an everywhere positive lapse and 
ends at spatial infinity in both universes. 

The above arguments do not mean that it is impossible 
to construct quasi-cquilibrium data with a combination 
of the conformal thin-sandwich approach and the punc- 
ture method. One might resign to constructing data in 
a slicing that cannot be everywhere stationary. For ex- 
ample, a construction might yield slicings that are quasi- 
stationary exterior to a neighborhood of the throat but 
become dynamical near the throat. The issue then is es- 
tablishing whether the resulting data are in fact in quasi- 
cquilibrium. Not enough is known yet to rule out the suc- 
cess of this technique. Alternatively one might consider 
allowing the lapse to cross zero. It is not yet completely 
clear what regularity conditions can be imposed or how 
the extrinsic curvature behaves in the neighborhood of 
these surfaces. At minimum this undertaking potentially 
involves a difficult numerical implementation at the crit- 
ical surfaces; the need to do so defeating the purpose of 
the puncture method. We conclude that the use of a 
combination of the conformal thin-sandwich and punc- 
ture methods for constructing quasi-equilibrium black 
hole initial data is less straightforward than it first ap- 
peared and may in fact be impractical. 
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